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DESIGN EVALUATION OF DISTRIBUTED SATURATION 
CHARACTERISTICS IN INTEGRATED DEVICES 

By W. W. Happ 

National Aeronautics and Space Administration 
Cambridge, Mass, 
and 

L. B. Dickson 
Motorola Semiconductor Co. 
Scottsdale, Arizona 


SUMMARY 

To evaluate the characteristics of an integrated circuit 
transistor operating in the saturated mode, it is necessary 
to determine the parameters for the distributed model. In 
addition, the lateral currents in the base and collector regions 
cause corresponding variations in the junction potentials. 
Formulas are given for the series collector resistance, and it 
is shown that normalized "A" parameters can be measured and will 
be constant for a given impurity profile. The difference equa- 
tions are solved with the aid of a computer. The solutions give 
not only the terminal saturation characteristics but also the 
internal current densities. Comparison of experimental values 
of the saturation characteristics with theory shows excellent 
agreement . 


INTRODUCTION 

A previous investigation (refs. 1,2) presented new 
approaches to model distributed saturation resistance in 
integrated devices and to develop a method for calculating 
the saturation voltage, specifically: 

(1) The transistor is separated into five distinct 
regions . 

(2) A distributed model is developed for each region 
which includes the effect of the base resistance and 
collector resistance. 

(3) A method for accurately calculating the value of 
these resistances using conformal mapping is given. 

This investigation will examine the validity and limita- 
tion of these modelling procedures, specifically: 

(1) The model parameters are determined from the geometry 
of the integrated device. 



(2) A specific integrated circuit transistor is analyzed 
numerically. 

(3) Experimental results are compared, with computer- 
aided design calculations. 

Finally, the use of these techniques and models in the 
design of integrated circuits is discussed in the light of 
current manufacturing experience and current laboratory 
prototypes now under development. 


DETERMINATION OF THE MODEL PARAMETERS 


One of the requirements of a useful transistor design 
model is that the relationship between the parameters and the 
geometry of the transistor being modeled can be determined by 
calculation or measurement. The model developed above has this 
desirable property. (Symbols used are defined in the Appendix.) 

Expressions for the A(e,c) parameters for drift transis- 
tors have been developed by Narud and Myer (ref. 3) and 
Callahan (ref. 4) . The dc expression for A^^ 
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where K accounts for the field and is the ratio of the base 
majority carriers at the emitter junction to the base majority 
carriers at the collector junction. Ag is the emitter area, 
and the other parameters have the standard notation (ref. 1) . 

The expressions for the other A(e,c) parameters and 
current gains are similar and involve the same terms (ref. 3). 

It is difficult to calculate the A(e,c) parameters with 
accuracy, mostly because of the unknown lifetime involved 
in Lp. However", the A(e,c) parameters per unit area may be 
accurately measured. Such measurements must be made at low 
current densities where the effects of the lateral current 
flow are negligible. 

Fortunately, it is a common practice for semiconductor 
manufacturers to use relatively few different impurity pro- 
files. An impurity profile defines the collector material 
and base and emitter diffusions. Generally, each manufacturer 
uses one standard inpurity profile to fabricate most saturat- 
ing logic circuits. If the A(e,c) parameters (normalized 
with respect to area) are measured for this standard impurity 
profile, the necessary A(e,c) parameter information is avail- 
albe and the model may be used to predict V^e(sat) 
various transistor geometries. 
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Figure 2.- Series collector 
resistance Rsc 
a function of 
geometry m/L and 
D/L from 

Ul _ tan h irL/4d 

U2 tan h (L/4d + TTm/2d) 


Approximation A>>1 in Table I is accurate to within 1% 
and leads to the simple equation for the calculation of the 
series collector resistance: 

R” = la ' fa >> 1 <2) 

This is the case when the distance between the emitter 
and collector contact is much greater than the thickness of 
the N region of the semiconductor material. Approximation 
4 in Table I is important when the widths of the emitter and 
collector contacts become small, and the contact resistance 
becomes quite large. Rearranging the results of- this approxi- 
mation, 

R ^ + .88 + £n fj (3) 

where 1/f = 1 - exp (— irm/2d) 

The last term in Eq. (3) becomes very large for small 
contacts, thus the equation predicts the effect of the con- 
struction or contact resistance. Equations (2) and (3) , then, 
give a straightforward method for calculating the series 
collector resistance in terms of the device geometry for the 
two specific cases. 


THE BASE RESISTANCE, Rg 

The sheet resistance (average resistivity divided by the 
width of layer) for thin diffused layers of silicon has been 
calculated by Irwin (ref. 6) , and his results are widely used 
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by most semiconductor manufacturers. Using these results, 
the sheet resistance (Psb^ base in Region I may be 

calculated. Then, the value of Rg for each section in Region 
I is simply: 




*^sb 


m/2 

Z 


where m is emitter width, n is the number of section, and Z is 
the emitter length. The base resistance in Region II, Rg , 
can be calculated using the formulas developed for R ^ 

and image techniques. An average resistivity, p, must be 
used. Its value is found by multiplying the total base 
sheet resistance by the base-collector junction depth. 


In this section, the methods for finding A(e,c) para- 
meters, Rg, Rgi/ and Rsq have been discussed. With this in- 
formation available, the model may be used to calculate the 
Vce(sat) integrated circuit transistor. 


A NUMERICAL EXAMPLE 


The saturation characteristics of the transistor, with 
dimensions given in Figure 3, will be evaluated from theory 
and compared with corresponding experimental results. 

Table II gives a summary of the measured A(e,c) parameters, 
alphas, and the calculated resistance. 



L = .40 mils 

L,= .30 mils 
b 

Z = 3 mils 

Figure 3.- 


C=Collector contact 
B=Ohmic base contact 
E=Emitter contact 


Scale drawing of 
the integrated cir- 
cuit transistor 


5 



TABLE II.- NUMERICAL VALUES OF PARAMETERS 
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9.3 S 
3.1 X 
.96 

.85 intrinsic - Region I 
14 ohms total 

.07 ohms/section in Region I 

.49 ohms/section in external regions 

4.17 ohms/section in Region I 

1.75 ohms/section in external regions 


From these values in Table II, the following can be 
observed: 

1) Although Age is larger than Aqq, the effective Aqq 
is less than the effective A^g, because of area 
considerations . 

The effective values are: 

A = 3 X 10~^^/mil^ 
ee ' 

-17 2 

A = 7 X 10 /mil 
cc ^ 

2) The emitter base foirward bias potential, Vg, is 
larger than the collector forward bias potential, 

V , at the same current level. 

The resistors are given in ohms per section and the 
intrinsic Region I of the model in Figure 4 will be divided 
into 20 sections. The value of Rgg for each section is 
found by taking the ratio of the length of one section under 
the emitter to the total length between the emitter and 
collector contacts and multiplying this ratio times Rg total. 
This procedure essentially linearizes Rgc appears to be 

in conflict with its method of calculation. However, as will 
be shown later, the results from this approximation are fairly 
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K = 0 K = I K = N 

Figure 4.- Distributed model of region I - 
the intrinsic region 

accurate. Dividing the extrinsic regions into 20 sections 
gives considerably larger values per section than in Region I, 
due to the larger sections. The same procedures have been 
followed in the base regions. 

COMPUTER EVALUATION - STRATEGY 

To gain insight into the operation of the integrated 
circuit transistor in saturation, the base and collector 
currents, as well as corresponding voltages, are calculated 
for each of the forty sections of the model. A high-current 
(IC = 40 ma) example and a low-current (I^ = 3 ma) example 
were examined . ^ 

Collector Voltage and Injected Current, JC 

Figure 5 shows the forward biased collector voltage V^, 
as a function of the lateral distance from the center of 
the transistor. For convenience, the cross-section of the 
transistor is reproduced at the top of the figure. To fully 
understand the results of Figure 5, one should re-examine the 
dimensions of the transistor given in the top view (Figure 3) . 
In the high -current mode, the differences in slope are evident, 
and the collector base junction is more heavily forward biased 
underneath the emitter and less foirward biased at the edge 
of the base than is the same junction for the low current 
case. The importance of this gradient in the collector bias 
is best appreciated if one realizes that for a change of 18 
millivolts, the current density, J^,, injected across the 
junction is reduced by one half. Then, the current being in- 
jected by the collector diode at X = .05 is more than two 
orders of magnitude greater than the current through the 
diode of X = . 40 for the high-current case. 
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X (mils) 

Figure 5.- Computer evaluation of the collector base 

forward bias potential as a function of distance 

The plot of Jcf shown in Figure 6, illustrates this fact. 
Clearly, the assumption that very little current is flowing 
in the base to the right of the base contact is a valid 
assumption for the high-current case. In the low-c\arrent 
case, the lateral current effects are less promounced, but 
the same trends hold. It is interesting to note that the 
diode at X = .40 is injecting more heavily for the low-current 
case than the high-current case. This is a result which 
would not be predicted by the classical models. 

Lateral Collector and Base Currents 

The total collector current (1^) flowing laterally in 
each section of the transistor (Figure 7) increases rapidly 
in the intrinsic region underneath the emitter and diminishes 
in the extrinsic base regions. The total base current flow- 
ing in each section of the model (Figure 8) also increases 
rapidly in Region I underneath the emitter, but continues to 
increase in the extrinsic regions. In fact, most of the base 
current originates in the extrinsic regions, and hence, comes 
from the collector. 
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Figure 6.- Computer eval- 
uation of the 
current injected 
from the collector 
into the base as 
a function of 
distance 


Figure 7.- Computer eval- 
uation of the 
total collector 
current flowing 
laterally as a 
function of 
distance 



Emmiter Current Density, Jg and Vjjg 

Hausser (ref. 7) has shown that Je is a function of the 
lateral distance (X) for the intrinsic region in the active 
mode. Figure 9 shows the variation of Jg for an integrated 
circuit transistor in the saturated mode. These results were 
obtained from the model consisting of 20 sections for the 
intrinsic region. The crowding effect or the change of cur- 
rent density with lateral distance is due to the potential 
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Figure 8.- Computer eval- 
uation of the 
total base current 
as a function of 
distance 


Figure 9 . - Computer evaluation 
of the injected 
emitter current as a 
function of distance 
for region I 



X (mils) 


drop of the lateral base current flow, and is much more 
pronounced at high currents. Figure 10 gives the variation 
of the base to emitter voltage in the intrinsic region. When 
the Ik total X Rbi extrinsic potential change is added to the 
junction voltage of the last section shown in the figure, the 
total potential between emitter and base contacts is foxand. 

No general curves of versus Ig were made, but the results 

of Figure 10 give a 10-mv agreement with the experimental 
results . 
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Figure 10.- Computer evaluation of the emitter 

base potential as a function of distance 

Effective 

All of the above results deal with currents and voltages 
inside the transistor rather than the terminal properties; 
however, as previously stated, these results may be used to 
calculate effective. The results of the calculations for 
the transistors at both current levels are given in Table III 

TABLE III.- COMPARISON OF CALCULATED AND EXPERIMENTAL VALUES. 


PARAMETER 

I TOTAL 
c 

(3.0 ma ) 

I TOTAL 
c 

(40 ma) 

"bi 


mamp 

.60 

3.67 

I . 

Cl 


mamp 

.73 

5.75 

«I 


calculated 

0.140 

0.360 



experimental 

0.148 

0.405 


% 


error 


6 % 


12 % 


SIGNIFICANCE OF RESULTS 


From the examination of these results and Figures 5 and 
6 , the following significant conclusions can be drawn: 

(1) At the high current levels, most of the carrier 
injection from the collector occurs near the emitter, 
thus causing a higher alpha inverse. 

(2) Differences in the injection level are due to the 
lateral pinch off of the collector base voltage as 
shown in Figure 5 . 

(3) The lateral pinch off at high -current levels causes 
an increase in the current injected from the collec- 
tor near the edge of the emitter, and thus, the 
assumption that only carriers injected in the in- 
trinsic region will be collected becomes less valid. 

(4) Because of the results of C, “leff will be somewhat 
larger at high current densities than the model 
predicts . 


^ce (sat) 

The experimental values of several transistors of 
Vce(sat) given by the curve of Figure 11. For collector 

currents of 3 milliamps, the model predicts the saturation 
voltage within 5 millivolts, and for collector currents of 
40 ma, the agreement is within 25 millivolts with the excep- 
tion of the measurements for low base currents. Both results 
are surprisingly accurate, but more important the general 
shape of the curves is in agreement with theory. As the base 
current is increased with a constant terminal collector 
current, the collector base junction becomes more heavily 
forward biased increases and decreases. 

^sc 

The series collector resistance as measured on the 
transistors was 8 ohms. This is the resistance from the two 
halves of the transistor in parallel and corresponds to a 
total collector resistance of 16 ohms for the proposed model, 
computed to a calculated value of 14 ohms. The error is 
probably due to the fact that RgQ was linearized, plus the 
fact that there is a finite (approximately .5 ohms) resistance 
in the connections to the integrated circuit transistor. 


12 



=39ma 




''ce(sat) 'b 

o = EXPERIMENTAL 
• = THEORETICAL 


1 ^ _i_ 

200 300 400 

(mv) 


Figure 11.- Comparison of computer evaluation with exper- 
imentally determined saturation characteristics 


CONCLUSIONS 


(1) The proposed model predicts the series collector 
resistance, the effective alpha inverse, and the saturation 
voltage for widely separated ' collector current levels, and 
just as important, it predicts the various current densities 
and voltages inside the integrated circuit transistor. 

(2) Parameters governing the operation of this transis- 
tor model were specified; and it was shown that, although 
these parameters are difficult to calculate, a set of normal- 
ized numerical values may be determined for a given impurity 
profile used to fabricate a type of integrated circuit. This 
is an important and favorable result, because semiconductor 
manufacturers tend to use only a few types of standard im- 
purity profiles to fabricate their circuits and each of these 
can be characterized. 


(3) An illustrative example has been presented for which 
the proper parameters have been calculated and used in the 
model. 
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APPENDIX 


A 


a 


n 


a 


I 


D 


I 

e 

J 


k 

kT 

q 


K,K^^ 


K 


K 

L 

1 

M 

N 


N 


1 


P 

no 


Glossary of Terms 

A basic parameter used to model transistors - defined 
in Eqs. (1) and (2) 

The common base current gain of a transistor operating 
in the normal mode, i.e. , otn = Ic/Ie, = constant 

The common base current gain of a transistor with the 
collector acting as the emitter, i.e., in the inverse 
mode 

Permitivity 

Diffusion coefficient of minority carrier in a semi- 
conductor 

Base current 

Collector current 

Emitter current 

2 

Current density in amps/cm 

Boltzman's constant = 8.63 X 10~^ eu/°K 

Reference voltage at 25 °C of value - 26 millivolts 

The complete elliptic integral of the first kind and 
its conjugate 

The field factor in Eq. (12) . K is the ratio of 
the doping in the base at the emitter to that at the 
collector 

The general section of the distributed mode 1 

Diffusion length of holes in N-type semiconductors 

The section immediately under the base contact 

The section at the edge of the emitter 

The section at the edge of the base-collector 
junction 

The equilibrium concentration of minority carriers in 
the base of a PNP transistor 
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q 


Electronic charge - 1.6 X 10 


-19 




Distributed base resistance under the emitter (in the 
active region) . 


% 


1 


Distributed base resistance in the external region. 


Rr 


1 


The lumped resistance external to the collector-base 
junction. 


R The total series collector resistance 

sc 

p Resistivity 

V .*The saturation voltage of an ideal one dimensional 

ce(sat) transistor with no series resistance 


V , The saturation voltage 
ce(sat) 

Emitter-base potential at any point 

V Forward biased collector-base potential at any point. 

^ The base width of the transistor in the active region. 
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